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Abstract 

The van der Waals (VDW) equation of state is a simple and popular model to describe the pressure 
function in equilibrium systems of particles with both repulsive and attractive interactions. This 
equation predicts an existence of a first-order liquid-gas phase transition and contains a critical point. 
Two steps to extend the VDW equation and make it appropriate for new physical applications are 
carried out in this paper: 1) the grand canonical ensemble formulation; 2) the inclusion of the quantum 
statistics. The VDW equation with Fermi statistics is then applied to a description of the system of 
interacting nucleons. The VDW parameters a and b are fixed to reproduce the properties of nuclear 
matter at saturation density no = 0.16 fm~ 3 and zero temperature. The model predicts a location 
of the critical point for the symmetric nuclear matter at temperature T c = 19.7 MeV and nucleon 
number density n c = 0.07 fm~ 3 . 
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I. INTRODUCTION 


The van der Waals equation of state is a simple analytical model of the pressure function p 
for equilibrium systems of particles with both attractive and repulsive interactions. The VDW 
model contains the first-order liquid-gas phase transition which ends at the critical point. In 
the canonical ensemble (CE), where independent variables are temperature T , volume V, and 
number of particles N, the VDW equation of state has most simple and transparent form, (see, 


e.g., Refs. 


> 2 ]), 


p(T, n) = 


NT 


N 2 


n T 




a n 


( 1 ) 


V - bN V 2 1 -bn 
where a > 0 and b > 0 are the VDW parameters which describe attractive and repulsive 
interactions, respectively, and n = N/V is the particle number density. The first term in the 
right-hand-side of Eq. (JT]) corresponds to the excluded volume (EV) correction, which manifests 
itself in a substitution of a total volume V by the available volume, V aw = V — bN. The second 
term comes from the mean field which describes attractive interactions between particles. In 
order to apply the VDW equation of state to systems with variable number of particles it is 
necessary to switch to the grand canonical ensemble (GCE). This procedure was first performed 
for the EV model, i.e., for a = 0 in Eq. (JTJ), in Refs. [3, 4]. In our recent paper p], the full 
VDW equation (JT]) , with both attractive and repulsive terms, was transformed from the CE to 
the GCE for systems with Boltzmann statistics. There are several physical situations when the 
GCE formulation is desirable (see Ref. j|j] for details). Note that the EV and VDW models can 
also be conveniently treated within the GCE in a framework of the thermodynamic mean-field 
approach (see Refs, [fj j|). 

Equation (JTJ) is valid for classical systems, where the effects of quantum statistics are ne¬ 
glected. In the present paper we suggest a generalization of the VDW equation to include effects 
of the quantum statistics. Proper treatment of quantum effects appears to be crucially impor¬ 
tant for a description of statistical equilibrium at small temperatures. The quantum statistics 
formulation is much easier to introduce in the GCE than in the CE. This is an additional 
physical example where the GCE formulation is particularly helpful. Thus, we use our recent 
results of the GCE formulation [fj] as a starting point for a quantum generalization of the VDW 
equation of state. As a next step, the VDW equation of state with Fermi statistics is used to 


2 










describe nuclear matter. The VDW parameters, a and b, which correspond, respectively, to 
attractive and repulsive interactions between nucleons, are fixed to reproduce the properties of 
the symmetric nuclear matter at zero temperature: saturation density Hq = 0.16 fm -3 , binding 
energy per nucleon, —16 MeV; and zero pressure, p = 0. 

The paper is organized as follows. In Sec. QT] the VDW equation of state is transformed into 
the GCE, and the quantum statistical formulation of this equation is elaborated in Sec. IIII1 In 
Sec. |IV] the VDW equation of state with Fermi statistics is applied to a description of nuclear 
matter. A summary in Sec. |V] closes the article. 


II. VDW EQUATION FOR THE BOLTZMANN STATISTICS IN THE GCE 


The VDW pressure function (JT]) corresponds to the Boltzmann approximation, i.e., the effects 
of quantum statistics (Bose or Fermi) are neglected. In our recent paper the VDW equation 
of state was formulated in the GCE. The GCE pressure, p(T, /i), is a function of temperature T 
and chemical potential p. It contains a complete information about thermodynamical functions 
of the system. Particle number density n(T, p), entropy density s(T, /a), and energy density 
e{T, p) can be presented in terms of p and its T and p derivatives: 






+J* 


dT 


dp 


P 


( 2 ) 


/' ' A 4 \~r~/ T 

For a = b = 0 the above VDW equations are reduced to the ideal gas expressions for classical 
particles. 

The VDW equation of state in the GCE is obtained in the form of a transcendental equation 
for particle number density n = n(T, p) as a function of T and p [5]: 

?r ld (T, p*) „ m bn 


n(T,p) = 


p* = p - T- 


1 + bn id (T,p*) ’ r ~ 1 — bn 

where n ld is a particle number density in the ideal Boltzmann gas, 


+ 2 an 


( 3 ) 


?r ld (T, p) = exp 


p\ dm 2 T 


2n 2 


Ko 


TJ ’ 


( 4 ) 


with d being the degeneracy factor and m the particle mass, K 2 (x) the Bessel function. Note 
that the relativistic form of a dispersion relation is considered, uo(k) = \Jm 2 + k 2 , where u> and 
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k are the free particle energy and momentum, respectively. The GCE VDW pressure p(T, //) 
is then obtained by inserting n(T,n) (J3J) into Eq. dTJ) . 

The VDW pressure (JTJ) is a unique function of variables T and n for all T > 0 and 0 < n < 
1/6. The VDW equation of state contains a first-order liquid-gas phase transition and has a 
critical point. The critical point (T c ,n c ) corresponds to the temperature and particle number 
density, where 


dp 

dn 


= 0, 


d 2 p 
dn 2 


= 0 


The thermodynamical quantities at the critical point are equal to 




_ 8a_ 
c 27 b ’ 


nc= 36’ Pc = 


27 b 2 


(5) 

( 6 ) 


At T > T c the following equation is always valid, 

dp 
dn 


> 0, 


while at T <T C the unstable region appears with 


dp 

dn 


< 0 . 


(7) 


( 8 ) 


T 


Therefore, the VDW isotherm p(n, T) at T < T c has a local maximum at n — ri\ and a local 
minimum at n = n 2 > n\. The unstable part (J8]) of the VDW isotherm at the interval [ 711 , 712 ], 
together with two additional parts - [n g ,ni\ and [ 712 , n-i] (they are called metastable) - are 
transformed to a mixture of two phases: a gas with density n g < ri\ and a liquid with density 
7 ii > n 2 ■ This is done according to the Maxwell rule of the equal areas (see, e.g., Refs. 2]), 
which leads to a constant pressure p(T, n g ) = p(T,ni ) inside the density interval [n g ,ni\. 

In the GCE the mixed phase region appears in a different way. At T > T c there is a unique 
solution of Eq. (]3|) . while at T <T C it may have either one solution or three different solutions 
for particle number density n(T,/j). Therefore, either one or three different solutions may also 
appear for the VDW pressure p(T,p). In a case when three different values of p(T,/i) are 
possible, the solution with a largest pressure survives in accordance to the Gibbs criterion (see 
Appendix [A] for details). The gas-liquid mixed phase in the T-/j plane belongs to the line 
p = p(T), where the solutions n g (T,p ) and ni(T,p ) correspond to equal pressures, p g (T,p ) = 


4 











The classical Boltzmann statistics leads to nonphysical behavior in zero temperature limit. 
This is already seen on an ideal gas level. For the ideal Boltzmann gas an entropy density in 
non-relativistic limit T/m 1 is equal to 


id c±L 
^Boltz 


n 


id 


T 


m 5r r 

m + —1 — 
2 



( 9 ) 


In Eq. (j9]) we use expressions p ld = n ld T and e ld = n ld {in + 3T/2) for the ideal gas pressure 
and (non-relativistic) energy density, respectively. Using an asymptotic expansion for the K 2 
Bessel function at large arguments, K 2 (x) = \Jnf{2x) exp(—x), one finds from Eq. (HD that to 
have a finite (nonzero) value n 0 of particle number density at T —y 0 the chemical potential 
should be equal to 

3T, /m/ . 2 tt n . . 

P = m -x- hi {T c 0 ) , c 0 = -— . (10) 

2 m 

Thus, only one limiting value, p = m, is admitted in the Boltzmann gas at T = 0 (this 
corresponds to zero value of chemical potential, /i n 0 n-rei = p — m, used in non-relativistic 
statistical physics). For p > m or p < m at T = 0 one finds for the particle number density 
n — 0 or n — oo, respectively. 

Therefore, the entropy of the ideal Boltzmann gas (19D at T —> 0 is 


s ld = n 0 


\ + ^ln(T/c 0 ) 


( 11 ) 


and it becomes negative in zero temperature limit, in a contradiction with the 3rd law of 
thermodynamics. The quantum statistics is needed to describe a physical system at T —y 0. 


III. VDW EQUATION OF STATE WITH QUANTUM STATISTICS 

Quantum generalization of the VDW equation of state is not a trivial task. Let us outline 
some general requirements for the quantum version of this equation of state: 

1. It should be transformed to the ideal quantum gas at a = 0 and 6 = 0. 

2. It should be equivalent to the classical VDW equation of state (HD in a region of thermo¬ 
dynamical parameters where quantum statistics can be neglected. 
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3. The entropy should be a non-negative quantity and go to zero at T —> 0. 


The pressure of the ideal quantum gas in the GCE reads 


p w (i» = i 


d 3 k 


k 2 


(2 71- ) 3 y/m 2 + k 2 


exp 


\/ m 2 + k 2 — p 


T 


+ rj 


-l 


( 12 ) 


In Eq. ()T2]h q equals +1 for Fermi statistics, -1 for Bose statistics, and 0 for the Boltzmann 
approximation. All other thermodynamical functions can be calculated from Eqs. (J2]). Ideal 
quantum gas expressions for thermodynamical functions satisfy the 3rd law of thermodynamics, 
i.e., s > 0 and s —* 0 at T —* 0. 

Let us now formulate a generalization of the VDW equation of state which includes effects 
of the quantum statistics. Note that p(T , p) for the Boltzmann case can be rewritten using 
Eq. ([3]) as 

p{T,fj) = p id (T, p*) — an 2 , (13) 


where 

H* — n — bp(T, p) — abn 2 + 2 an . (14) 


The function p ld in Eq. (TT3|) corresponds to the ideal gas pressure in the Boltzmann approxi¬ 
mation, i.e. rj = 0 in Eq. (1T2]) . We suggest the quantum VDW equation of state in the same 
form as Eq. (1T3]) but with ideal quantum gas pressure p ld , i.e., for quantum case we propose to 
take rj — ±1 in Eq. f| 12 p . which corresponds to the Fermi or Bose statistics. 

In accordance with ([2]), one has the following for the particle number density: 

= (IOt = " ld(T ’ , ‘' )(1 " H ( 1 + 2 “ti) -2 “"ir (15) 

This equation can be transformed to 


n id (T, p,*) (1 - bn) - n 


dn\ 

1 + 2a— = 0 . 


(16) 


The solution of this equation, which has a physical meaning, reads 
n id (T, p*) 


n(T,p) = 


1 + bn id (T, p*) 


n ld (T,p) = 


d 

2iT 2 


dkk- 


exp 


V m 2 + k 2 — p 


+ p 


-l 


. (U) 


and it has the same form as Eq. (J2]). However, a principal difference is that n ld in Eq. ((171) is 
a particle number density of the ideal quantum gas, whereas n ld in Eq. (J3J) corresponds to the 
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ideal classical gas, i.e., r/ = 0, and is given by Eq. We also note that, in the quantum case, 
expression (fT4l) for the shifted chemical potential p* should be used instead of ([3]). 

Equations (U3l) and (1X71) correspond to the system of two equations for two unknown func¬ 
tions: p(T, p) and n(T, p). The VDW model defined by these equations possesses all the 
required properties. First, at a = 0 and b = 0 Eqs. (flTSD and (1171) are reduced to the ideal 
quantum gas expressions. Second, for those T and /i* values, where quantum expressions for p ld 
and n ld can be approximated by the Boltzmann statistics, i.e., by Eq. © for n ld and p ld = Tn ld 
for the ideal gas pressure, Eqs. m and 0 EP become automatically equivalent to the classical 
VDW equation of state (0Q). Third, the entropy density has the following form 


s(T, p) 


dp\ _ s id (T, p*) 

- 1 + bn id (T, p*) ’ 


(18) 


thus, it is always positive for the quantum, ideal gas expressions of s ld , and s —* 0 at T —> 0. 
The energy density can be calculated from Eq. 02]) as 


e{T,p) 


e id (T, p*) 


an 


n , 


(19) 


where e id (T, p) is the average energy per particle in the ideal gas, 

w* = ' (20) 

One can rewrite the VDW pressure as a function of temperature T and particle density n. 
It follows from Eq. (Eli that p* can be written as a function of T and n, 

^ • < 21 > 

where /i ld (n, T) is the chemical potential of the ideal quantum gas, which is a solution of the 
following transcendental equation for the given n and T : 

n = 2^1 dkk 
Equation (fT3l) can be then rewritten as 


exp 


\fm 2 + k 2 — p 

Y~ 


id 


n -1 


+ v 


( 22 ) 


P 



n 

1 — bn 



— an 2 . 


(23) 
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One can easily check that this equation coincides with Eq. (JT]) in a case of the Boltzmann 
statistics and, thus, can be indeed regarded as a quantum generalization of the classical VDW 
equation in the CE. It can be also instructive to consider a formulation of the VDW equation 


with quantum statistics within the thermodynamic mean-held approach developed in Refs. 
8 ]. This is presented in Appendix lESl 


IV. NUCLEAR MATTER 


In this section the VDW equation of state with quantum statistics is used to describe the 
properties of symmetric nuclear matter. Namely, a Fermi gas of nucleons (m = 938 MeV and 
d = 4) is considered with attractive and repulsive interactions described by the a and b VDW 
parameters, respectively. A study of nuclear matter has a long history. The thermodynamics 
of nuclear matter and its applications to the production of nuclear fragments in heavy ion 


13] in 1980s. A review of these early developments can be 


collisions were considered in Refs. 

found in Ref. [mJ]. Nowadays, the properties of nuclear matter are described by many different 


models, particularly by those which employ a self-consistent mean held approach .MM 
Excluded-volume corrections in the mean-held models have been considered in Refs. 4, 7, 20]. 
Experiment ally, a presence of the liquid-gas phase transition in nuclear matter was hrst reported 


in Refs. 


21 


23] by indirect observations. The hrst direct measurements of the nuclear ca’ 


oric 


curve were done by the ALADIN collaboration 24], later followed by other experiments 


25 


26]. 


Our consideration will be restricted to small temperatures, T < 30 MeV, thus, a pion 
production will be neglected. In the present work, we also neglect a possible formation of 
nucleon clusters (i.e., ordinary nuclei) and baryonic resonances (like N* and A), which may 
be important at low and high baryonic density, respectively. Within these approximations, the 
number of nucleons N becomes a conserved number and an independent variable in the CE. 
The chemical potential /i of the GCE regulates the number density of nucleons. 






























A. Properties at T = 0 


For calculations of the thermodynamic functions in the GCE, Eqs. (1T31) . (fT7]h (fT8]h and (TT9l) 
will be used. In terms of variable p* (TT5|1 . thermodynamical functions of the quantum VDW 
gas can be presented in terms of the corresponding functions of the ideal quantum gas as the 
following: 


£(T,/i) = 


n' d (T,p*) 

1 + bn id (T, p*) 

e id (T,p*) 

1 + bn id (T, p*) 


= p ld (T, p*) - a 


— a 


n ld (T, p*) 


n lA (T,p*) 

1 + bn id (T, p*) 

s id (T, p*) 


s(T,n) = 


1 + b n iA (T, p*) 


1 + b n iA (T, p*)_ 

At zero temperature the ideal gas quantities in Eqs. (1241) and (1251) can be written as 


n ld (T = 0, p*) = 


p ld (T = 0, p*) = 


d 

2^ 


d 

6tt 2 


dkk 2 = (/i* 2 - m 2 ) 3/2 , 

07 T z 




dk 


d 


VW+ 


rrr 


487T 2 


p* a/ p* 2 — m 2 (2/i* 2 — 5m 2 ) 


167T 2 


rrr In ■ 


m 


A 


+ 


^2 _ ^^2 


\/ V* 2 -r 


e id (T = 0, p*) = 


d 


dkk 2 Vk 2 + m 2 — “ - p* Wp* 2 — m 2 (2 p* 2 — m 2 ) 

lo7r- 


s Id (T = 0, p*) = lim 
t-> o 


d 

2^ 

d m 4 i n . _ ) 

167T 2 n* 4- yjfi *2 — m 2 ’ 

e id (T, p*) + p id (T, p*) - p* n id (T, p*) 


(24) 

(25) 


(26) 


(27) 


m 


T 


= 0 . 


(28) 

(29) 


We fix parameters a and b in such a way to reproduce the properties of nuclear matter in its 
ground state (see, e.g., Ref. 27]), i.e., it should be p = 0 and £/n = m + Eb = 922 MeV at 
T — 0 and n — no = 0.16 fm -3 . Here Eb = —16 MeV is the binding energy per nucleon. One 
then finds, a = 329 MeV fm 3 and b = 3.42 fm 3 . Note that parameter b of the proper particle 
volume can be expressed in terms of hard-core radius r as b — 167rr 3 /3. This gives r = 0.59 fm 
for the hard-core nucleon radius. 
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Figure 1: Dependence of pressure p (a) and binding energy Eb (b) on the nucleon density n at 
T = 0. The VDW parameters are a = 329 MeV-fm 3 and b = 3.42 fm 3 (r = 0.59 fm). The open 
circle corresponds to the ground state of nuclear matter. The dash-dotted line corresponds to the 
metastable part of VDW isotherm, whereas the dotted line corresponds to the unstable part. 




Figure 2: The same as in Fig. [I] but for small values of n. 

The pressure p and binding energy Eb as functions on nucleon density n at T — 0 are shown 
in Fig. m (a) and (b), respectively. The stable VDW isotherms are depicted in Fig. ID by solid 
lines, while the metastable and unstable parts are depicted by dash-dotted and dotted lines, 
respectively. At very small densities a gaseous phase with almost ideal gas behavior is always 
present. At T = 0 this phase, seen more clearly in Fig. [2J can, however, exist as a metastable 
state only. 

Note that at any T > 0 the chemical potential has a well-defined limiting behavior /i —y —oo 
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at n —> 0. At T = 0 the situation is different: at n = 0 the chemical potential may have any 
value smaller than the particle mass. The mixed gas-liquid phase at T = 0 is depicted by the 
horizontal lines in Fig. |T] (a) and (b). The two coexisting phases at T = 0 are the liquid phase 
with n = 0.16 fm -3 , p = 0, and /i = 922 MeV, and the gaseous phase with n — 0, p — 0 
and ji = 922 MeV. This corresponds to the Gibbs conditions of phase equilibrium, i.e., equal 
temperatures, pressures, and chemical potentials for coexisting phases. The stable gaseous 
phase at T = 0 is, in fact, a vacuum with n = 0. 


B. Phase diagram 


The VDW pressure isotherms are depicted in (T, v) and ( T,n ) coordinates ( v = 1 /n) in 
Fig. [3] (a) and (b), respectively. They are calculated within the quantum VDW equation of 
state using Eq. (12cS|) with a = 329 MeV fm 3 and b = 3.42 fm 3 . The critical temperature is 
found to be T c = 19.7 MeV. The value of the critical temperature in our model is close to the 
experimental estimates in Refs. [25, [26|. At T < T c two phases appear: the gas and liquid 
phases separated by a first-order phase transition. The mixed phase region is obtained from 
the Maxwell construction of equal areas for p{y) isotherms (see Appendix [Aj) , and it is depicted 
by horizontal lines in Fig. [3] (a) and by the shaded grey area in Fig. [3] (b). The nucleon number 
density at the critical point is found to be n c = 0.07 fm -3 = 0.4 no- Normal nuclear matter 
with n = no = 0.16 fm -3 and T = 0 corresponds to a point placed exactly on the boundary 
between the mixed and liquid phases. Note also that the maximal value of the nucleon number 
density in the VDW model is n max = 1/6, which is equal to n max = 0.29 fm -3 for the chosen 
value of parameter b. 

In the mixed phase region, the particle number density is given by 


n = Cn g + (1 - 0 n i > ( 30 ) 

where £ and 1 — £ are the volume fractions of the gaseous and liquid components, respectively. 
The values of n g and ni in Eq. (I5H are the particle densities of, respectively, the gaseous and 
liquid phases at the corresponding boundaries with the mixed phase. A behavior of the mixed 
phase at T = 0 is rather special. The stable gaseous phase is absent at T = 0, i.e., n g = 0 at 


11 







Figure 3: Pressure isotherms in (a) (p, v ) and (b) (p, n) coordinates, calculated in the quantum Van 
der Waals equation of state with parameters a = 329 MeV-fm 3 and b = 3.42 fm 3 (r = 0.59 frn). The 
dashed-dotted lines present the metastable parts of the VDW isotherms at T < T c , whereas the dotted 
lines correspond to unstable parts. The full circle on the T = T c isotherm corresponds to the critical 
point, while the open circle at T = 0 in (b) shows the ground state of nuclear matter. Shaded grey 
area in (b) depicts the mixed phase region obtained from the Maxwell construction of equal areas for 
p(v) isotherms in (a). 

the mixed phase boundary as this boundary starts from the point T = 0 and n = 0. Therefore, 
only a metastable gaseous phase at small densities can exist at T = 0 as depicted in Fig. [2] The 
stable gaseous phase exists however at any T > 0 for small enough values of the particle number 
density, smaller than the n g density of gaseous phase in the mixed phase region resulting from 
the Maxwell construction. 

Parameters of the critical point found in the VDW case with Fermi statistics for nucleons 
differ significantly from those values for the classical VDW gas. With the same VDW param¬ 
eters a and b as in the Fermi statistics, the classical VDW equation ([I]), i.e., with Boltzmann 
statistics, would give T c = 8a/27b = 28.5 MeV and n c = 1/3 b = 0.10 fm -3 . This further 
indicates an importance of the effects of quantum statistics: these effects are not only crucial 
in the limit T —> 0 but also remain quantitatively important even near the critical point. 

In Fig. [4]the phase diagram of the symmetric nuclear matter in (T, p) coordinates is depicted. 
The nucleon density at different temperature T and chemical potential p is presented. At T <T C 
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there is the T-p region with three different solutions for p(T, p) at given T and p. According to 
the Gibbs criterion, a solution with a largest pressure survives (see Appendix 0, and only this 
solution for n(T,p) is depicted in Fig. |4j The location of the critical point is shown in Fig. [4] 
by the solid circle, while the normal nuclear matter state corresponds to the open circle. Note 
that /i 0 = 922 MeV corresponds to the chemical potential of the normal nuclear matter which 
is placed on a boundary with the liquid phase. The values p < po are forbidden at T = 0 (these 
values of p lead formally to n — 0). The values p > po at T = 0 are possible and correspond 
to the nuclear liquid. 

The phase transition line, p = /i m j x (T), shown in Fig. SI starts from the normal nuclear 
matter state with T = 0, p = 922 MeV and ends at the critical point with T c = 19.7 MeV, 
Pmix(T c ) = 908 MeV. This line presents the whole mixed phase region shown by the grey area 
in Fig. |3j At each T <T C two solutions, n g (T , p) and ni(T , p), with different particle densities, 
n g (T, p) < ni(T,p ), and equal pressures, p g (T } p) = pi(T,p), exist at the phase transition line 
p = p c {T). On this line, the discontinuities of thermodynamical quantities n, s, and s take 


place. 

At T > T c there is only one solution n(T, p) for any T and p values, i.e., there are no distinct 
gaseous or liquid phases. Nevertheless, as seen from Fig. [3] very rapid, although continuous, 
changes of particle number density take place in a narrow T-p region even at T > T c . This is 
a manifestation of the so-called smooth crossover phenomenon. 

At any T > 0, there are no restrictions on possible values of the chemical potential; i.e., any 
values of p between — oo and +oo are possible. When p decreases the particle number density 
decreases too and goes to zero at p —» — oo. At very small n, both the particle interactions and 
Fermi statistics effects become negligible. The system of nucleons behaves then as the ideal 
Boltzmann gas. In an opposite limit, p —> oo, nucleon density n(T, p) goes to its upper limiting 
value 1/6. The VDW pressure behaves then approximately as p = nT/( 1 — bn) and goes to 
infinity. Different theoretical models and their comparison with experimental estimates of the 
nuclear matter properties have been widely discussed in the literature (see, e.g., Refs. 128 — 
ft. In the present paper we do not attempt to make any detail comparison of the developed 
VDW quantum model with existing data for nuclear matter. Some extensions of the model will 
probably be needed. These questions are, however, beyond the scope of the present paper. 
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Figure 4: Density of the symmetric nuclear matter in (T, //) coordinates, calculated in the quantum 
van der Waals equation of state with parameters a = 329 MeV fm 3 and b = 3.42 fm 3 (r = 0.59 fm). 
The open circle denotes the ground state of nuclear matter, the full circle corresponds to the critical 
endpoint, and the phase transition curve is depicted by the solid black line. 

V. SUMMARY 

In the present paper we have formulated a generalization of the van der Waals equation of 
state to include effects of quantum statistics. In the grand canonical ensemble a system of two 
transcendental equations for the pressure and particle density is obtained. These equations can 
be solved for all possible values of temperature, T > 0, and chemical potential, —oo < /i < oo. 
Our quantum generalization of the VDW equation satisfies all basic requirements: it reduces 
to the ideal Fermi or Bose gas for a = 0 and b = 0, to the classical VDW equation in the 
Boltzmann limit, and it satisfies the 3rd law of thermodynamics, i.e. s —> 0 as T —> 0. 

The VDW equation with Fermi statistics has then been applied to a system of interacting 
nucleons to describe the properties of symmetric nuclear matter. The VDW parameters a and 
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b of interacting nucleons are fixed by the properties of the nuclear matter ground state: T — 0, 
p — 0, n — no = 0.16 fm -3 , and Eb = —16 MeV. We find a = 329 MeVfrn 3 and b = 3.42 fm 3 . 
With these parameters the VDW model predicts a first-order liquid-gas phase transition with 
a critical endpoint located at T c = 19.7 MeV and n c = 0.07 fm -3 . Extensions of the presented 
formulation as well as new physical applications will be the subject of further studies. 

Acknowledgments 

We would like to thank M. Gazdzicki, A. G. Magner, S. Mrowczynski, and K. Rcdlich for 
fruitful comments and discussions. This work was partially supported by HIC for FAIR within 
the LOEWE program of the State of Hesse and by the Program of Fundamental Research of 
the Department of Physics and Astronomy of National Academy of Sciences of Ukraine. 


Appendix A: Gibbs criteria and Maxwell construction 


A statistical system in the GCE is defined by two independent variables, T and p. Two 
distinct phases, gas and liquid, coexist if their pressures are equal, p g (T,p ) = pi(T,p). In the 
case of p g (T,fi ) ^ pi(T,p ), only a phase with a larger pressure survives. These statements are 
known as the Gibbs criteria for the the first-order phase transition (see, e.g., Refs. Q,3). We 
now prove that the Gibbs criteria are equivalent to the Maxwell construction of equal areas for 
the VDW equation of state. 

At T < T c , the Maxwell construction replaces a part of the VDW isotherm p = p(v, T) 
by the horizontal line p = p mix , which corresponds to the mixed phase region for all v in the 
interval [vi,v g ]. This is shown in Fig. 0 The Maxwell equal areas are 


rvo 


dv \p D 


p(v,T )] = / dv \p(v,T) 


Pn 


(A.l) 


fvi 


'VQ 


where Vi < v 0 < v g and p(v 0 ,T ) = p mix . The replaced parts of the isotherm are interpreted as 
metastable (dp/dv < 0) and unstable (dp/dv > 0) states. They are shown in Fig. [5] by the 
dashed-dotted and dotted lines, respectively. 
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Using the thermodynamical identity 



one can present the chemical potential (ia at any point A on the isotherm as 


(A.2) 


/u 


CPA 


Hb 


dp' v(p', T) 


'Pb 


(A.3) 


where B is an arbitrary point on the isotherm, and integration in (IA.3P is performed along the 
path from point B to point A on the isotherm. With Eqs. (I A. 1(1 and (I A. 3(1 one can easily prove 



Figure 5: A generic VDW isotherm p(v, T ) for T < T c . Points l and g on the isotherm correspond to 
the boundaries of the mixed phase. In (a) point 1 (2) on the isotherm corresponds to the maximum 
(minimum) value of the chemical potential inside the mixed phase, which is reached at the boundary 
between metastable and unstable phases, while point V (2') corresponds to the same value of the 
chemical potential reached in the pure liquid (gaseous) phase. In (b) points a\ (gaseous phase), 02 
(metastable liquid phase) and 03 (unstable phase) which all correspond to the same value of the 
chemical potential are depicted on the VDW isotherm. 

that 

p(vi,T) = n{v g ,T ) = // mix , (A.4) 

thus, the Maxwell and Gibbs constructions for the mixed phase are equivalent. 

According to Eq. (IA.3D the chemical potential decreases with v if ( dp/dv)T < 0, and increases 
if (dp/dv)T > 0. Therefore, inside the mixed phase region [vi,v g ] the chemical potential reaches 
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its minimal value p — p 2 < p mix at v = v 2 and its maximal value p — p 1 > p m ix at v — v \. The 
points v\ and v 2 , where ( dp/dv)x = 0, correspond to the boundaries between the metastable 
and unstable parts of the VDW isotherm. At v < Vi and v > v g the chemical potential is 
a monotonously decreasing function of v, with p — > 00 at v —> b and /1 —> —00 at v —> 00 . 
Therefore, there is a point v 2 > v g in the gaseous phase where p = /i 2 , and a point v[ < Vi in 
the liquid phase where p = p\. These two points are depicted in Fig. [5] (a). At both /i > /i\ 
and ji < n 2 the GCE VDW pressure p(T, / 1 ) is a unique function. On the other hand, there are 
three different solutions for the VDW pressure at /X 2 < A 4 < Mi- 

Let us first consider AD < ha < hmix- There are three points on the VDW isotherm with 
H = /ia shown in Fig. 0 (b): point cq in the gaseous phase with v g < < v 2 , point a 2 in 

the metastable liquid phase with with Vi < v 2 < u 2 , and point a 3 in the unstable phase with 
V‘2 < v 3 < Vi . Using Eq. (IA.3[) one Ends that 

PPniix. PPniix. /*Pmix PP3 

/ dp'v(p',T ) = / dp'v(p',T) = / dp'v(p',T) — / dp'v(p',T) . (A.5) 

dpi dpi d P 2 Jp 2 

Applying the mean value theorem to integrals in Eq. (IA.5[) one obtains 

(Pmix - Pi) Vl = (pmix - pV) V 2 = (pmix ~ P 2 ) E 31 - (p 2 ~ P%) T 32 , (A.6) 

where Vi <v 2 < v 2 , v g <Vi< Vi < v 3l < v 2 < v 32 < v 3 , and p 2 = p(v 2 , T). It follows from 
Eq. (1A.6P that > p 2 and p\ > p%] i.e., at p 2 < p < /i rn ix the gaseous phase should be realized 
according to the Gibbs criterium as its pressure is larger than pressures of both metastable and 
unstable states. 

The same arguments are applied to A^mix < p < pi and show that the liquid pressure is 
then larger than the pressures of both metastable and unstable states. Therefore, for the 
VDW equation of state the Maxwell construction of the equal areas and Gibbs criteria are fully 
equivalent. The Maxwell construction is applied in the CE, whereas the Gibbs criteria are used 
in the GCE. We emphasize that this statement is valid not only for the classical VDW equation 
of state, but also for the VDW equation with Fermi statistics with isotherms depicted in Fig. [3j 
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Appendix B: Thermodynamic mean-field approach 

In the framework of the thermodynamic mean-field (TMF) approach [§-8] the pressure and 
particle number density are presented as 

p(T,fx) = p id [T,p id (n,T)] + P ex (n, T) , (B.l) 

n{T,n) = n id [T, p — U(n, T)] , (B.2) 

where P ex (n, T) and U(n,T) are, respectively, the excess pressure and the thermodynamic 
mean field. The presence of nonzero quantities P ex (n, T) and U(n,T) in Eqs. (IB. ID and (IB. 21) 
correspond to interaction between particles, and the condition of thermodynamic consistency 
reads 8]: 

dU dP ex . 

n ~dn = ~BrT ' (B 3) 

For a specific choice of P ex (n, T ) and U(n, T) functions, one proceeds by solving of Eq. (1B.2[) for 
n = n(T, u), and then the pressure p(T, p) can be obtained from Eq. (IB. II) (see some examples 
in Ref. |8(). 

The VDW equation of state with quantum statistics defined by Eqs. (fT31) and (1171) can be 
rewritten in the TMF form (IB.II) and (IB. 2D . In order to determine P ex (n, T) and U(n,T) we 
rewrite Eq. (123|) for the VDW pressure as 


P = p ld p> ld (n,T)] + p ld [T,p 
Comparing (IB.ID and (IB.dD one finds 


.id 


n 


1 — bn 


, T 


- p ld [I> M (n, T)} - a 


.id 


n 


pex 

-M/dW 


(' n,T ) = p ld T,p 


.id 


n 


,T 


- p‘ d [T,^(n,T)] - a 


. 1 — bn 

The mean field Uv$w(n,T) can be then calculated from Eq. (IB. 3D as 


n 


Pvdw(n,T) = 


1 dP vdw( n ',T) 


n 


dn' 


dn' . 


(B.4) 


(B.5) 


(B.6) 


For the Boltzmann statistics, Eqs. (1B.5D and (IB.6D are simplified to the following analytical 
expressions: 

bn 0 . 

— an 2 , (B.7) 


Pyd W ( n , T) = Tn 


U, 


VdW 


[n,T) = T 


1 — bn 
bn 


1 — bn 


Tln(l — bn) — 2 an. 


(B.8) 
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